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ABSTRACT 

Qnasi-periodic oscillations (QPOs) at freqnencies near 1000 Hz were recently dis¬ 
covered in several x-ray binaries containing nentron stars. Two sonrces show no cor¬ 
relation between QPO freqnency and sonrce connt rate (Berger et ah 1996, Zhang et 
ah 1996). We snggest that the QPO freqnency is determined by the Keplerian orbital 
freqnency near the marginally stable orbit predicted by general relativity in strong grav¬ 
itational helds (Mnchotrzeb-Czerny 1986, Paczyhski 1987, Klnzniak et ah 1990). The 
QPO frequencies observed from 4U 1636-536 imply that the mass of the neutron star is 
2.02 ± 0.12 Mq. Interpretation of the 4.1 keV absorption line observed from 4U 1636- 
536 (Waki et ah 1984) as due to Fe xxv ions then implies a neutron star radius of 
9.6 ± 0.6 km. 


Subject headings: accretion, accretion disks — gravitation — relativity — stars: indi¬ 
vidual (4U 1636-536, 4U 1608-52) — stars: neutron — X-rays: stars 
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1. Introduction 

Observations of the low-mass x-ray binaries 4U 1608- 
52 (Berger et al. 1996) and 4U 1636-536 (Zhang et 
al. 1996) with the Rossi X-Ray Timing Explorer 
(Bradt, Rothschild, & Swank 1993) have revealed 
quasi-periodic oscillations (QPOs) at 800-900 Hz which 
show no or little dependence of QPO frequency on 
source count rate. In 4U 1608-52, the QPO frequency 
changes by less than 10% as the count rate changes 
by a factor of 5 (Berger et al. 1996). This is in 
sharp contrast to the strong correlation of QPO fre¬ 
quency and source count rate seen in the low-mass 
x-ray binaries 4U 1728-34 (Strohmayer et al. 1996) 
and 4U 0614-k091 (Ford et al. 1997); in 4U 0614-k091, 
the QPO frequency changes by more than a factor of 
2 as the count rate changes by a factor of 2 (Ford et al. 
1997). In 4U 1608-52 and 4U 1636-536, the QPO fre¬ 
quency is, at most, only weakly influenced by the lu¬ 
minosity or mass accretion rate of the source. We sug¬ 
gest that the QPO frequency is determined by the Ke- 
plerian orbital frequency near the marginally stable 
orbit predicted by general relativity in strong gravi¬ 
tational fields (Muchotrzeb-Czerny 1986; Paczyhski 
1987; Kluzniak & Wagoner 1985; Kluzniak et al. 
1990). 

2. Accretion Flow at the Marginally Stable 
Orbit 

In general relativity, no stable orbits exist close to a 
sufficiently compact massive object. The radius of the 
marginally stable orbit r^s is determined by the mass 
and rotation of the star and is independent of the 
mass accretion rate or radiation flux in the accretion 
disk. In a Schwarzschild spacetime, rms = QGM/c^, 
where M is the mass of the star. For ‘soft’ nuclear 
matter equations of state, the marginally stable orbit 
can lie outside the neutron star (Kluzniak & Wagoner 
1985; Cook, Shapiro, & Teukolsky 1994). If the neu¬ 
tron star’s magnetic field is sufficiently small, then 
the accretion disk will pass through r^s- Keplerian 
orbits near have frequencies near 1000 Hz and 
could generate QPOs which are independent of source 
count rate. 

An analysis of thin accretion disks by Muchotrzeb- 
Czerny (Muchotrzeb 1983; Muchotrzeb-Czerny 1986) 
showed that a stationary flow can not be supported 
near rms for viscosities greater than a critical value. 
Stationary flow breaks down at the sonic radius, 
where the radial velocity of matter in the accretion 


disk changes from subsonic to supersonic. This insta¬ 
bility may lead to observable oscillations either via 
Doppler beaming or eclipses and provides a mecha¬ 
nism to generate QPOs at Keplerian orbital frequen¬ 
cies near (Paczyhski 1987; Kluzniak et al. 1990; 
Miller, Lamb, & Psaltis 1996). 

The sonic radius is located close to rms, but may 
move due to variations in surface density in the in¬ 
ner disk (Muchotrzeb-Czerny 1986). The analysis 
of Muchotrzeb-Czerny (1986) suggests that the sonic 
radius will oscillate over a relatively narrow range 
near rms- If the surface density variations are ir¬ 
regular, then changes in the QPO frequency will be 
irregular. The QPO frequencies in 4U 1608-52 and 
4U 1636-536 show an irregular time evolution, wan¬ 
dering over a small range in frequency (50 Hz) on a 
time scale of 1000 s (Berger et al. 1996; Zhang et 
al. 1996). The observed range of QPO frequency 
variations imply variations in the orbital radius of 
Ar/r = 2Ai'/ii/ ~ 0.03. This is consistent with the 
allowed range for sonic radius variations calculated by 
Muchotrzeb-Czerny (1986) for a ~ 0.2 — 0.3, where 
a is a dimensionless viscosity parameter (Shakura & 
Sunyaev 1973). 

Motion of the sonic radius should occur on the 
viscous time scale t,, ~ R^/v where R is the disk 
radius and v is the kinematic viscosity. Using an 
a-prescription for the viscosity, we find that ti, ~ 
where Ai is the Mach number of the flow 
defined as the ratio of the circular velocity to the 
sound speed, and is the orbital time scale (Pringle 
1981). If the QPOs are associated with Keplerian or¬ 
bits then ~ 10“^ s. For the inner edge of the disk 
at a temperature kT ~ 1 keV, we estimate A4 ~ 300. 
For this model to hold, we must have a viscosity 
greater than the critical value discussed above, which 
imples a > Oc ~ 0.03. Estimates of a for the in¬ 
ner accretion disk are in the range 0.1 — 1 (Eardley 
& Lightman 1975). Viscosities in this range lead to 
viscous time scales of 10^ — 10^ s. This matches the 
observed 10^ s time scale for QPO frequency varia¬ 
tions in 4U 1608-52 and 4U 1636-536. 

We conclude that motion of the sonic point near 
the marginally stable orbit produces small QPO fre¬ 
quency variations which are consistent in time scale 
and frequency range with those observed in 4U 1608- 
52 and 4U 1636-536. The count rate independence 
and approximate constancy of the QPO frequency 
arises naturally from the count rate independence of 
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In one observation of 4U 1636-536, two QPOs are 
present at frequencies ofll71±llHz and 898.7 ± 
1.4 Hz (van der Klis et al. 1996). This is similar to 
the high-frequency QPOs in 4U 1728-34 (Strohmayer 
et al. 1996) and 4U 0614-1-091 (Ford et al. 1996) 
which have been interpreted with models in which 
the higher frequency is the Keplerian frequency in the 
inner region of the accretion disk while the lower fre¬ 
quency ‘beat’ signal results from differential rotation 
of the inner disk and the spinning neutron star (Al- 
par & Shaham 1985; Lamb et al. 1985). We suggest 
that the QPOs near 800-900 Hz in 4U 1608-52 and 
4U 1636-536 are beat frequency signals. Comparison 
with 4U 1728-34 and 4U 0614-091 supports this iden¬ 
tification. In these two sources the fractional RMS 
amplitude of the higher frequency QPO decreases 
with count rate while the amplitude of the lower fre¬ 
quency peak increases. Extrapolation of these trends 
would lead one to anticipate that only the beat fre¬ 
quency QPO should be visible in higher luminosity 
sources such as 4U 1608-52 and 4U 1636-536. 

A possible explanation of why only the beat fre¬ 
quency signal is visible is that the instability at the 
sonic point produces a modulation in the beaming 
pattern of the radiation but not in the total lumi¬ 
nosity of the source, while a luminosity variation at 
the beat frequency can result from radiation feedback 
from the neutron star onto the disk (Miller, Lamb, & 
Psaltis 1996). If the optical depth near the neutron 
star is large, then the observed amplitude of beam¬ 
ing modulations will be reduced relative to luminosity 
modulations (Kylafis & Phinney 1989), and, for high 
optical depths, only the beat-frequency QPO would 
be visible. If the optical depth between the star and 
the disk is sufficiently large, the radiation feedback, 
and therefore the beat-frequency QPO, may be sup¬ 
pressed (Miller, Lamb, & Psaltis 1996). 

The strong dependence of QPO frequency in 4U 1728- 
34 and 4U 0614-091 on source intensity implies that 
the accretion disks in these sources must be disrupted 
before reaching r^s- The qualitiatively different be¬ 
havior of the QPOs in 4U 1728-34 and 4U 0614-091 
versus those in 4U 1636-536 and 4U 1608-52 is a nat¬ 
ural consequence of where the inner disk is disrupted 
relative to r^s- 

3. Discussion 

If, indeed, the QPOs in 4U 1608-52 and 4U 1636- 
536 are associated with Keplerian orbits near the 


marginally stable orbit, then they provide us with 
a tool for determination of the masses of the neu¬ 
tron stars in these systems (Kluzniak et al. 1990). 
For 4U 1636-536, the beat frequency of 899 Hz, ob¬ 
served simultaneously with the Keplerian frequency 
of 1171 Hz, lies near the maximum beat frequency 
observed. Since the analysis of Muchotrzeb-Czerny 
(1986) indicates that the sonic radius is always equal 
to or larger than rms for viscosities above the critical 
value, this suggests that the corresponding Keplerian 
orbital frequency of 1171 Hz corresponds to an or¬ 
bit near r^ns- The analysis of Chen and Taam (1993) 
also shows that the sonic radius lies very close to, al¬ 
though perhaps inside, rms for accretion rates below 
the Eddington rate. 

In a Kerr spacetime, the relation between the mass 
of the star and the Keplerian orbital frequency is M = 
2.198MQ{iyK/W00 Hz)-i(l - 0.748j)”\ to first order 
in j, where j = clto/GM'^ is a dimensionless mea¬ 
sure of the angular momentum of the star (Bardeen, 
Press, & Teukolsky 1972; Kluzniak et al. 1990). The 
stellar rotation frequency a;/27r = 272Hz, and an as¬ 
sumed moment of inertia / = 2.0 x lO^^gcm^, imply 
a neutron star mass of 2.02Mq. Allowing an error of 
50 Hz in the Keplerian frequency, sufficient for both 
the experimental and theoretical uncertainties, con¬ 
tributes an error of 5% to the mass estimate. Allow¬ 
ing moments of inertia in the range 1 — 3 x 10^^ g cm^ 
contributes an an uncertainty of 4% to the mass esti¬ 
mate. We conclude that the mass of the neutron star 
in 4U 1636-536 is 2.O2±O.12M0. This mass is within 
the allowed range of neutron star masses (Kalogera 
& Baym 1996) and is not unreasonable given that ac¬ 
cretion onto the neutron star has occured. 

It is interesting to note that an absorption feature 
at 4.1 ±0.1 keV has been observed in thermonuclear x- 
ray bursts from 4U 1636-536 (Waki et al. 1984). Sim¬ 
ilar absorption features have been seen from 4U 1608- 
52 (Nakamura, Inoue, & Tanaka 1988) and EXO 1747- 
214 (Magnier et al. 1989) by Tenma and EXOSAT. 
The most plausible interpretation of the line-like ab¬ 
sorption feature is absorption by Fe XXV ions in mat¬ 
ter accreted onto the neutron star surface after the 
outburst (Waki et al. 1984; Inoue 1988). The trans¬ 
verse Doppler effect is negligible for a neutron star 
spin rate of 272 Hz. In this case, the redshift of the 
line and the neutron star mass quoted above implies 
a neutron star radius of 9.6 ± 0.6 km. This result 
is consistent with the mass-radius relation for the 
AV14±UVH equation of state (Wiringa, Fiks, & Fab- 
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rocini 1988) and allowing for rotation (Cook, Shapiro, 
& Teukolsky 1994). We note this equation of state 
produces a moment of inertia of 2.0 x lO^^gcm^ for 
a 2.0Mq star, consistent with our assumed I above. 

If confirmed, the relation between high-frequency 
QPOs in x-ray binaries and the marginally stable or¬ 
bit in general relativity may open a new experimental 
venue for the study of strong gravitational fields and 
will permit the determination of neutron star masses 
as described above. Further work is required to im¬ 
prove our theoretical understanding of accretion flows 
near r^s and to obtain more detailed observations of 
high-frequency QPOs and the x-ray binary systems 
that produce QPOs. Independent determination of 
the mass of the neutron star in 4U 1636-536, via mea¬ 
surement of the binary orbital parameters, would be 
of great interest. Comparison of the mass with the 
radius of the marginally stable orbit determined from 
the QPO frequency would permit a test of strong field 
and inertial frame dragging effects in general relativ¬ 
ity (Kluzniak et al. 1990). 

We thank M.C. Miller, F.K. Lamb, and M. Ruder- 
man for useful discussions. 
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